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ABSTRACT 

We investigate the well-measured spectral and energetic properties of 20 
gamma-ray bursts (GRBs) in their cosmological rest frames. We find a tight 
relation between the isotropic-equivalent 7-ray energy E iso , the local peak energy 
E' of the vF v spectrum, and the local break time t' b of the GRB afterglow light 
curve, which reads E iso f b oc E'™ 5 * - * (xl = 1-40; fl M = 0.27, fi A = 0.73). Such 
a power-law relation can be understood via the high-energy radiation processes 
for the GRB prompt emission accompanying the beaming effects. We then con- 
sider this relation as an intrinsic one for the observed GRB sample, and obtain 
a constraint on the mass density Qm = 0.241q 12 (1°") f° r a ACDM universe, 
and a Xdof = 1-33 for Qm « 0.3 and fl\ ~ 0.7. Ongoing GRB observations in the 
Swift era are expected to confirm this relation and make its cosmological utility 
progress much. 

Subject headings: gamma rays: bursts — cosmology: observations — cosmology: 
distance scale 



1. Introduction 



Gamma-ray bursts (GRBs) are the most powerful explosions in the universe since the 
Big Bang. Their cosmological origins are identified by the redshift measurements of their 
exploded remnants, usually named "GRB afterglows", or their host galaxies. GRBs are 
common regarded as jetted phenomena, supported by the observational evidence that an 
achromatic break appears in the afterglow light curve, which declines more steeply than in 
the spherical model largely due to the edge effect and the laterally spreading effect (Rhoads 
1999; Sari et al. 1999) or the off-axis viewing effect (Rossi, Lazzati & Rees 2002; Zhang & 
Meszaros 2002a; Berger et al. 2003). With the advantages of huge energy release for the 
prompt emission and immunity to dust extinction for the 7-ray photons, GRBs are widely 
believed to be detectable out to a very high redshift of z ~ 10 — 20 (Lamb & Reichart 2000; 
Ciardi & Loeb 2000; Bromm & Loeb 2002; Gou et al. 2004). 
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Similar to the Phillips relation in Type la supernovae (SNe la; Phillips 1993), a tight 
relation in GRBs, linking a couple of energetic and observational properties, could make 
GRBs new standard candles. Amati et al. (2002) reported a power-law relation between 
the isotropic-equivalent 7-ray energy E iso of a GRB and its rest-frame vF v peak energy 
E' . Unfortunately, the large scatter around this relation stymies its cosmological purpose. 
However, after joining the burst's half-opening angle 9, Ghirlanda et al. (2004a) found that 
a tight relation is instead between the collimation-corrected energy and E' p , which reads 
Ej = £i so (l — cos#) oc E' p R (k is the power). Recently, Liang & Zhang (2005) reported a 
multi- variable relation between E iso , E' p) and the local achromatic break t' b of the burst's 
afterglow light curve, i.e., E iso oc E' p K1 t' b K2 (ki and k 2 are the powers). Applying these two 
relations to different observed GRB samples, meaningful cosmological constraints have been 
performed in a series of works (e.g. Dai, Liang & Xu 2004; Ghirlanda et al. 2004b; Firmani 
et al. 2005; Xu, Dai & Liang 2005; Ghisellini et al. 2005; Qin et al. 2005; Mortsell & 
Sollerman 2005; Liang & Zhang 2005; Bertolami & Silva 2005; Lamb et al. 2005 for a GRB 
mission plan to investigate the properties of dark energy). 

In this paper we investigate the well-measured spectral and energetic properties of the 
up-to-date 20 GRBs in their cosmological frames and report a new relation between Ex 
(E x = E iso t' b ) and E', which is Ex oc E' p (where b is the power). More importantly, 
inspiring constraints on cosmological parameters can be achieved if this newly found relation 
is indeed an intrinsic one for the observed GRB sample. 



2. Data Analysis and Statistical Result 

The "bolometric" isotropic-equivalent energy of a GRB is given by 

E iso = 4nd 2 L S^k(l + z)- 1 , (1) 

where S* 7 is the fluence in the observed bandpass and k is a multiplicative correction relating 
the observed bandpass with a standard rest-frame bandpass (1-10 4 keV) (Bloom, Frail & 
Sari 2001), with its fractional uncertainty 




Thus, the collimation-related energy E x (see §4 of this work) is 

E x = E iso t' b = A7rd 2 L S 7 kt b (l + z)~\ (3) 
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where % is the observed break time of the afterglow light curve in the optical band, with its 
fractional uncertainty 

(ifHf) 2 +(?) 2 +fe) 2 (4) 

To calculate /c-correction of a burst, one is required to know its spectral parameters 
fitted by the Band function, i.e., the low-energy spectral index a, the high-energy spectral 
index f3, and the observed peak energy E p = E'/(l + z) (Band et al. 1993). In this paper, 
we considered those GRBs with known redshifts, fluences, spectral parameters, and break 
times, and thus collected a sample of 20 bursts shown in Table 1. The criterions of our data 
selection are: (1) The fluence and spectral parameters of a burst are chosen from the same 
original literature as possible. If different fluences are reported in that literature, we choose 
the measurement in the widest energy band. If this criterion is unsatisfied, we choose the 
fluence measured in the widest energy band available in other literature; (2) The observed 
break time for each burst is taken from its data in the optical band; (3) When the fractional 
uncertainties of E p , and k are not reported in the literature, they are taken to be 20%, 
10% and 5%, respectively; (4) To be more reliable, a lower limit of the fractional uncertainty 
of t b is set to be 10%. 

In this section, we investigate the E iso — E' p and Ex — E' relations for a Friedmann- 
Roberston- Walker cosmology with mass density VLm and vacuum energy density Q\. There- 
fore, the theoretical luminosity distance in equations (1) and (3) is given by 

d L = c(l + ^)^ 1 jn fc |- 1 / 2 sinii{jn fc | 1 / 2 

x f dz[(l + z) 2 (l + n M z)-z(2 + z)n A ]' 1/2 }, (5) 
Jo 

where £l k — 1 — Vl M — ^a, Ho = 100 h kms _1 Mpc _1 is the present Hubble constant, and 
"sinn" is sinh for £l k > and sin for tt k < 0. For Q k = 0, equation (5) degenerates to be 
c(l + z)Hq 1 times the integral (Carroll et al. 1992). According to the conclusions of the 
Wilkinson Microwave Anisotropy Probe (WMAP) observations, we here choose Om = 0.27, 
Qa = 0.73, and h = 0.71 (e.g. Spergel et al. 2003). The dashed and solid lines in Figure 
1 respectively represent the best-fit powerlaws for the Ei SO — E' and Ex — E' relations 
(weighting for the uncertainties on both coordinates, see Press et al. 1999). We find 

E- 

10 52 ergs 
with a reduced x 2 — 5-99, and 

Ex 
10 52 ergs 



EL 



x 1.90±0.06 



10 0.27±0.03 V__ | (6) 

1 10 2 KeV/ 



I n-0.04±0.04 I E p . 

l 10 2 KeV j {) 



1.95±0.08 
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with a reduced x 2 — 1-40. As can be seen, although the power of the E iso — E' p relation in this 
work is roughly consistent with that in Amati et al. (2002) derived from 9 BeppoSAX bursts, 
i.e., E' p oc ^0 o 52±0 - 06 (-^2 _ o.91), the dispersion around this relation increases seriously. 
Instead, the scatter of the Ex — E' relation is rather small. One may notice that it is the 
Ex — E' p relation rather than the Amati relation whose power more converges at ~ 2. Based 
on the above statistical findings, we consider the Ex — E' p relation for the observed GRB 
sample to constrain cosmological parameters. 



3. Constraints on Cosmological Parameters 

The E x — E' p relation can be given by 

£ x /10 52 ergs =10 a (^/100KeV) fe , (8) 

where the dimensionless parameters a and b are assumed to have no covariance. Combining 
equations (3) and (8), we derive the apparent luminosity distance as 

d L = 9.142 x 10 a/2 (£;/100) 6/2 (l + z){kS y t b )~ 1 / 2 Mpc, (9) 

with its fractional uncertainty being 

( as A 2 1 ( ak \ 2 1 ( at A 2 1 ( b ° E * 

d L ) ~~ \2Sj \2k) \2t b ) \2 E p 

+ (!^(!^)\ 

where all the quantities are assumed to be independent of each other and their uncertainties 
follow Gaussian distributions. The distance modulus is calculated by \i = 5 log(<iL/Mpc)+25, 
and its uncertainty is computed through = (5/ In 10) {crd L /d>L)- 

Although the Ex — E' relation is similar to the Phillips relation in SNe la, the methods 
by which to constrain cosmological parameters are different (Riess et al. 1998; Firmani et al. 
2005, Xu et al. 2005). For SNe la, a Phillips-like relation is known to those relatively high-z 
objects after it has been calibrated by nearby well-observed SNe la, because the theoretical 
luminosity distance is irrelevant with cosmological parameters, i.e. = z(c/H ), when 
z 1. While for GRBs, one won't know that sole set of (a, b)- parameters in the E x — E' p 
relation until a \ow-z GRB sample is established 1 . Therefore, the x 2 statistic for GRBs is 



20 



x 2 (n,a,b\h) = J2 



k=l 



fith(z k ; tt\h) - Ho hs (z k ; fl, a, b\h) 

°MobsOfc; ^ a, b, a a /a, a b /b) 



2 



(11) 



Possible cosmic evolution for this relation herein cautioned. 
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where Q = (Qm, ^a) denotes a certain cosmology, and h is taken as 0.71 in this work. 

We carry out a Bayesian approach to obtain GRBs' constraints on cosmological param- 
eters. For clarity, we divide it into three stages with the detailed steps as follows: 

Stage I 

(1) fix Qi = (Q M ,Q\)i, (2) calculate /i th and E x for each burst for that cosmology, (3) 
best fit the Ex — E' relation to yield (a, b)i and (cr a /a, <?b/b)i, (4) apply the best-fit relation 
parametrized by (a,b)i and (a a /a, 0&/&)j to the observed sample, and thus derive /z b s and 
o" Mobs for each burst for that cosmology, (5) repeat Steps 1—4 from i — 1 to i — N to obtain 
A*th, /^obs and <T Mobs for each burst for each cosmology; 

Stage II 

(6) re-fix Qj, (7) calculate x 2 {^j\^i) by comparing // t h(%) with // bs(^i), 0> O bs(^)> 
and then convert it to a conditional probability, i.e., probability for Qj which is contributed 
by the relation calibrated for by the formula of P(Qj\Qi) oc exp[— % 2 (fi J -|Jl i )/2], (8) 
repeat Step 7 from i = 1 to % = N to obtain an iterative probability for cosmology Qj by 
P lte (Qj) oc J^exp[— x 2 (r2j|r2j)/2] x P im (fij) (here the initial probability for each cosmology 

i 

is regarded as equal, i.e., P im (Q) = 1), (9) repeat Steps 6—8 from j — 1 to j — N to obtain 
an iterative probability P ltc (Q) for each cosmology; 

Stage III 

The iterative probability P lte (f2) for each cosmology derived on Step 9 is no longer 
equal to its initial probability P mi (fi), but it has not reached the final/converged probability 
P fin (fi). So the following process is to (10) assign P itc (fi) on Step 9 to P ini (fi) on Step 8, then 
repeat Steps 8—9, and thus reach another set of iterative probabilities for each cosmology, 
(11) run the above iteration cycle again and again until the probability for each cosmology 
converges, i.e., P ltc (f2) =^ P fin (f2) after tens of cycles. 

In this method, to calculate the probability for a favored cosmology, we consider contri- 
butions of all the possible Ex — E' p relations associated with their weights. The conditional 
probability P(ll, |f2j) denotes the contribution of some certain relation, and P fin (f2j) weights 
the likelihood of this relation for its corresponding cosmology. Therefore, this Bayesian 
approach can be formulized by 

N / N 

P fin (^) = ^P(^I^) x P fin (fi t ) / ]TP fin (^) (j = 1,N). 
i=i / i=i 

Shown in Figure 2 are the constraints for (fi M ,fi A ) from 20 observed GRBs (solid 
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contours), using the likelihood procedure presented here. The GRB dataset is consistent 
with Q M « 0.3 and ~ 0.7, yielding a Xdof = 1-33. For a flat universe, we measure 
tt M = 0.24+°;}^ at the 68.3% confidence level (C.L.). The best-fit point is Q M = 0.24 and 
Q\ = 0.76 (red cross). Shape of the elliptic confidence contours trends to well constrain Qm 
and provide evidence for cosmic acceleration with an enlarged GRB sample. For comparison, 
we also plot the constraints from 157 gold SNe la in Riess et al. (2004) (dashed contours). 
Closed tt M — ranges derived from these two datasets at the same C.L.s have conveyed 
the advantages of high- ,2 distance indicators in constraining cosmological parameters as well 
as other cosmological issues. Additionally, GRBs are complementary to SNe la on the 
cosmological aspect. We see from the figure that a combination of SNe and GRBs wound 
make the cosmic model of Qm ~ 0.3 and Q A ~ 0.7 more favored, thus better in agreement 
with the result form WMAP observations (e.g. Bennett et al. 2003). GRBs are hopeful to 
be new standard candles. 



4. Conclusions and Discussion 

We report a new relation between GRB rest-frame energetics and spectra, which is 
E iso t' b /10 52 ergs = 10- - 04±0 - 04 (E;/100KeV) L95±0 - 08 ( X l = 1-40) for Q M = 0.27, fi A = 0.73 
and h = 0.71. Considering this power-law relation for the 20 observed GRBs, we find a 
constraint on the mass density Q M = 0.24+°;}^ (la) for a flat universe together with a 
XL = 1-33 for n M ~ 0.3 and Q A « 0.7. 

As previously discussed, the Ghirlanda relation reads E 1 = E iso (l — cos#) oc E' p K (k is 
the power) under the framework of uniform jet model. Within this model, one can calculate 
the half-opening ang le of a GRB jet by 9 oc ^{n^^Erl 1 ^ where n is the uniform 
circumburst medium density, and r? 7 denotes the conversion efficiency of the initial ejecta's 
kinetic energy to 7-ray energy release (Rhoads 1999; Sari et al. 1999). Thanks to strong 
collimation for GRBs, i.e. 9 <C 1, the Ghirlanda relation becomes E 1 oc (E iso ^) 3 / 4 (n^ 7 ) 1 / 4 . 
In Ghirlanda et al. (2004a), the term of to? 7 was assumed to be highly clustered for the 
whole sample of 15 bursts presented there, thus the Ghirlanda relation turns out to be E 1 oc 
{E iso t' b ) 3 / A . So if E iso t' b oc E' p b (b ~ 2), there will be E 1 oc E' p K (k ~ 1.5). In this sense, the 
Ey — E' and Ex~E' p relations are consistent with each other. However, 7] 1 should be different 
from burst to burst (e.g. < 1% — 90%), and n is variable for a burst in the wind environment 
and may expand in several orders (see Friedman & Bloom 2004 and references therein). 
Additionally, these two observables and their uncertainties are very difficult to be reliably 
estimated. Therefore, if the Ghirlanda relation is used to measure cosmology, the resulting 
constraints would depend on different input assumptions of r] 7 and n. For contrast, such 
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difficulties have been circumvented when applying the E x — E' relation. Moreover, Liang & 
Zhang (2005) reported a generalized relation between GRB rest-frame spectra and energetics 
could be written as E iso /I0 52 ergs = (0.85 ± 0.21) x (E' p /100 KeV) L94±ai7 x (t' b /l day)" L24±a23 
for a flat universe of Om = 0.28, using a multiple variable regression method for the 15 bursts 
presented there. Their conclusion is consistent with our statistical result (see §2 in this 
work). More loose constraints, however, were performed when the uncertainties of the fitted 
parameters in this relation were included into the error of the apparent distance modulus (see 
Fig 11 in Liang & Zhang 2005). Among the discussed three relations, the E x — E' p relation 
is the simplest and the latter two have the advantage of making the relations explicitly 
model-independent and eliminating the needs to marginalize over the unknown r/ 7 and n. 

So what is the underlying theoretical basis for the E x —E' p relation? At the present stage, 
plausible explanations mainly include: the standard synchrotron mechanisms in relativistic 
shocks (e.g. Zhang & Meszaros 2002b; Dai & Lu 2002), the high-energy emission form off- 
axis relativistic jets (e.g. Yamazaki et al. 2004; Eichler & Levinson 2004; Levinson & Eichler 
2005), and the dissipative photosphere model producing a relativistic outflow (e.g. Rees & 
Meszaros 2005). The scaling relation depends on the details of each model, and it resembles 
the observational result under certain simplification. Although models are different, the 
relation they support has made GRBs towards more and more standardized candles. 

XD is thankful to Z. G. Dai, E. W. Liang and B. Zhang for valuable discussions and D. 
Band in explaining the observations of GRB050525a. This work is supported by the National 
Natural Science Foundation of China (grants 10233010 and 10221001), and the Ministry of 
Science and Technology of China (NKBRSF G19990754). 
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Table 1: 

Sample of 20 7-ray bursts 



GRB 


Redshift 






S 7 K) b 


Bandpass b 


*bK) c 


References' 1 






KeV 




10~ 6 erg cm- 2 


KeV 


day 


(E p , S 7 , t h ) 


970828... 


0.9578 


297.7[59.5] 


-0.70, -2.07 


96.0[9.6] 


20-2000 


2.2(0.4) 


1,2,3 


980703... 


0.966 


254.0[50.8] 


-1.31, -2.40 


22.6[2.26] 


20-2000 


3.4(0.5) 


1,2,4 


990123... 


1.600 


780.8(61.9) 


-0.89, -2.45 


300.0(40.0) 


40-700 


2.04(0.46) 


5,5,6 


990510... 


1.619 


161.5(16.0) 


-1.23, -2.70 


19.0(2.0) 


40-700 


1.57(0.16) 


5,5,7 


990705... 


0.8424 


188.8(15.2) 


-1.05, -2.20 


75.0(8.0) 


40-700 


1.0(0.2) 


5,5,8 


990712... 


0.4331 


65.0(10.5) 


-1.88, -2.48 


6.5(0.3) 


40-700 


1.6(0.2) 


5,5,9 


991216... 


1.020 


317.3[63.4] 


-1.23, -2.18 


194.0[19.4] 


20-2000 


1.2(0.4) 


1,2,10 


011211... 


2.140 


59.2(7.6) 


-0.84, -2.30 


5.0[0.5] 


40-700 


1.56(0.16) 


11,2,12 


020124... 


3.200 


120.0(22.6) 


-1.10, -2.30 


6.8[0.68] 


30-400 


3.0(0.4) 


13,13,14 


020405... 


0.690 


192.5(53.8) 


0.00, -1.87 


74.0(0.7) 


15-2000 


1.67(0.52) 


15,15,15 


020813... 


1.255 


212.0(42.0) 


-1.05, -2.30 


102.0(10.2] 


30-400 


0.43(0.06) 


13,13,16 


021004... 


2.332 


79.8(30.0) 


-1.01, -2.30 


2.55(0.60) 


2 -400 


4.74(0.47) 


17,17,18 


021211... 


1.006 


46.8(5.5) 


-0.805,-2.37 


2.17(0.15) 


30-400 


1.4(0.5) 


19,19,20 


030226... 


1.986 


97.1(20.0) 


-0.89, -2.30 


5.61(0.65) 


2 -400 


1.04(0.12) 


17,17,21 


030328... 


1.520 


126.3(13.5) 


-1.14, -2.09 


36.95(1.40) 


2 -400 


0.8(0.1) 


17,17,22 


030329... 


0.1685 


67.9(2.2) 


-1.26, -2.28 


110.0(10.0) 


30-400 


0.48(0.05) 


23,23,24 


030429... 


2.658 


35.0(9.0) 


-1.12, -2.30 


0.854(0.14) 


2 -400 


1.77(1.0) 


17,17,25 


041006... 


0.7160 


63.4[12.7] 


-1.37, -2.30 


19.9(1.99] 


25-100 


0.16(0.04) 


26,26,27 


050408... 


1.2357 


19.93(4.0) 


-1.979,-2.30 


1.90(0.19] 


30-400 


0.28(0.17) 


28,28,29 


050525a. 


0.606 


78.8(4.0) 


-0.987,-8.839 


20.1(0.50) 


15-350 


0.20(0.10) 


30,30,30 



Note. — (a) The spectral parameters fitted by the Band function. The fractional uncertainty of E p is 
taken as 20% when not reported, and the fractional uncertainty of ^-correction is fixed as 5%; (b) The 
flucnccs and their errors in the observed energy band. The fractional uncertainty of S 7 is taken as 10% when 
not reported. The fluence and spectral parameters of a burst are chosen from the same original literature as 
possible. If this criterion is unsatisfied, we choose the fluence measured in the widest energy band available 
in other literature; (c) Afterglow break times and their errors in the optical band. A lower limit of 10% is 
set for the fractional uncertainties of 4; (d) References in order for E p [\a,0\), 5 7 (Bandpass), and t\,. 

References. — (1) Jimenez et al. 2001; (2) Bloom et al .2003; (3) Djorgovski et al. 2001; (4) Frail et al. 
2003; (5) Amati et al. 2002; (6) Kulkarni et al. 1999; (7) Stanek et al. 1999; (8) Masetti et al. 2000; (9) 
Bjornsson et al. 2001; (10) Halpcrn et al. 2000; (11) Amati 2003; (12) Jakobsson et al. 2003; (13) Barraud 
et al. 2003; (14) Berger et al. 2002; (15) Price et al. 2003a; (16) Barth et al. 2003; (17) Sakamoto et al. 
2004; (18) Holland et al. 2003; (19) Crew et al. 2003; (20) Holland et al. 2004; (21) Klose et al. 2004; (22) 
Andersen et al. 2003; (23) Vanderspek et al. 2004; (24) Price et al. 2003b; (25) Jakobsson et al. 2004; (26) 
Butler et al. 2005; (27) Stanek et al. 2005; (28) Sakamoto et al. 2005; (29) Godet et al. 2005; (30) Blustin 
et al. 2005. 
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Fig. 1. — The Ex/E iso — E' plane. Circles represent Ex vs. E' for 20 GRBs, while squares 
denote E- lso vs. E' p . Solid and dashed lines are the best fits of the Ex — E' and E iso — E' p 
relations, using Q M = 0.27, £7a = 0.73, and h = 0.71. 
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Fig. 2. — Joint confidence intervals (68.3%, 90%, 99%) in the Qm — plane from 20 
observed GRBs (solid contours), and from 157 gold SNe la (dashed contours). Red and 
black crosses mark the best fits of the two samples, while black dot denotes the concordance 
model of Qm = 0.3. 



